Does a low solar cycle minimum hint at a weak upcoming cycle? 
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ABSTRACT 

The maximum amplitude (R m ) of a solar cycle, in the term of mean sunspot numbers, is 
well-known to be positively correlated with the preceding minimum (i? mm )- So far as the long 
term trend is concerned, a low level of i? m i n tends to be followed by a weak R m , and vice versa. In 
this paper, we found that the evidence is insufficient to infer a very weak Cycle 24 from the very 
low Rmin in the preceding cycle. This is concluded by analyzing the correlation in the temporal 
variations of parameters for two successive cycles. 



Subject headings: sun: activity — sun: general — sunspots 



1. Introduction 

Studying the correlation between the maximum 
amplitude (i? m ) of a solar cycle and the preceding 
minimum (R m - m ) is useful for understanding the 
long-term evolution of solar activity. This can pro- 
vide information about the activity level of an en- 
suing cycle. The positive correla tion between R m 
and R m in is a well-known fact (jHathawav et al. 
20021) . As a natural consequence, the very low 



level of solar activity at the present time (around 
the onset o f Cycle 24) seems to be followed by a 
very weak (jSvalgaard et al.ll2005t ISchatterJ[2005h . 



or even the weakest cycle (|Lil 120091 ) . However, 
a lower R m i n has not always been followed by a 
weaker cycle. For example, a sm all R min precedes 
the g reatest R m in Cycle 19 ( Wang fc Sheelevl 
2009). Therefore, what information we can infer 
from the preceding minimum is worth re-analyzing 
carefully. We ask whether and how past cycles af- 
fect the present cycle. 

A more accurate prediction of solar activity 
is an important task in solar physics and space 
weather. Knowing in advance the activity level of 
an upcoming cycle is helpful in the launching and 
operation of spacecrafts. An underestimate of the 
activity level for the next cycle may let down our 
guard. One aim of this study serves to remind the 
space flight mission planners that they still need 



to remain vigilant to avoid unexpected troubles. 

In the present study, we use the 13- month run- 
ning mean of Zurich relative sunspot numbeiQ 
from January 1749 to April 2010 to determine 
the maximum (i? m ) and the preceding minimum 
(-Rmin) of the solar cycle. The correlation between 
R m and i? m i n for different periods of time is shown 
in Section [2j In Section |3l we examine the vary- 
ing trends of R m and R m i n using a quantity to de- 
scribe whether a parameter increases or decreases. 
Then, in Section 21 we analyze the temporal varia- 
tion in the correlation coefficient between R m and 
Rmin with a moving time window of five cycles. 
The results are briefly discussed and summarized 
in Section [5] 

2. Correlation between R m and i? m i n 

The parameters of R m and -R m i n since cycle n = 
1 are listed in Table [IJ and shown in Fig. la. 

It can be seen from Fig. la that R m and i? m i n 
have a similar long-term variation behavior: a 
lower (higher) level of i? m i n tends to be followed 
by a weaker (stronger) R m . Their correlation co- 
efficient is r = 0.56 at a confidence level (CL) of 
99%. The scatter plot of R m against i? m i n (tri- 



1 http: / / www.ngdc.noaa.gov/stp/SOLAR/ftpsunspot- 
numbcr.html 
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n(r= 0.18, CL=57%) 

Fig. 1. — (a) i? m (solid) and i? m i„ (dotted) since cycle 
n =1. (b) Scatter plot of R m against i? m i n (triangles), 
(c) Varying trend: V m (solid) and Vmin (dotted). 



angles) is shown in Fig. lb. Their least-squares-fit 
linear regression equation is: 



R m = 77.9 + 5.98-Rn 



a = 33.5, 



(1) 



where a is the standard deviation of the equation. 
Substituting the present value of i? m i n (1-7) into 
this equation, the peak of Cycle 24 is predicted as 
R m (24) = 88.0 ± 33.5 (labeled by a star), regard- 
ing the l-<7 uncertainty. When using the modern 
era data since Cycle 10, the peak of Cycle 24 is 
predicted to be higher, as i? m (24) = 98.4 ± 35.0. 
When using only the recent nine cycles since Cycle 
15, an even higher value of i? m (24) = 122.3 ± 36.5 
will be predicted for Cycle 24. 

However, [l| (|2009h inferred a rather low level 
for Cycle 24, i? m (24) = 58.0 ± 26.6, from a rela- 
tio nship of R m = 48 .8 + 5.39fl min ±26.6 derived 
by I Hathaway et all ( 2002 ). The data used to de- 
rive this relationship are those that are smoothed 
with the 24-month Gaussian filter, rather than 
the 'standard' 13-month running mean. So, the 
present minimum value (i? m i n = 1.7) of the 13- 
month running mean sunspot number is inappro- 
priate for use in inferring an i? m (24) value from 
this relationship. Besides, the i? m i n value in terms 
of the 24-month Gaussian filter is unknown within 
12 months of the minimum. (Even if this value is 
known, the result inferred from the above relation- 
ship has a different meaning.) 

We return to Fig. la. If using only the param- 
eters in the earlier cycles of n = 1-14 (left to the 



vertical line in Fig. la), the correlation coefficient 
between R m and i? m in increases to r(l-14) = 0.72 
at the 99% level of confidence. In contrast, for the 
recent cycles of n = 15-23, the correlation coeffi- 
cient is only r(15-23) = 0.23, which is statistically 
insignificant (CL < 50%). Therefore, the positive 
correlation between R m and -R m i n (0.56) is mostly 
contributed by the earlier cycles. The recent cy- 
cles, especially for cycles 15-19, seem to behave 
differently from the earlier cycles. It is then nec- 
essary to analyze whether the temporal variation 
in the correlation affects the future i? m value. 

3. Trends of variations in R m and i? m i n 

It should be noted in Fig. la that, for an indi- 
vidual cycle, the increase or decrease of R m does 
not always follow that of i? m i n . For example, i? m in 
decreases while i? m increases from n = 18 to 19. 
To demonstrate the behavior of increasing or de- 
creasing, we define the varying trends of i? m and 

-^min &S. 

V m (n) = Sgn(i? m (n) - R m (n - 1)), 
V min (n) = Sgn(i? min (n) - i? mi „(n - 1)), 

(2) 

where y = Sgn(x) is the sign function: y = 1 if x > 
0, y = -1 if x < and y = if x = 0. V m (n) = 
+1 refers to an increase in R m (n) > R nl (n — 1), 
and V m (n) = — 1 refers to a decrease in R m {n) < 
Rm(n — 1), and so on. V^n) = Vmi n (n) refers 
to the same trend of R m (n) and i? m i n (n). The 
values of V m (n) and V m in(n) are listed in Tabled] 
and shown in Fig. lc. For all cycles n = 2-23, 
there are 13/9 pairs of V m (n) and V m i n (ri) with the 
same/opposite trends. Their correlation is very 
weak, at r — 0.18, and statistically insignificant 
(CL = 57%). 

For the earlier cycles of n = 2-14, there 
are 10/3 pairs of V m (n) and V m - m (n) with the 
same/opposite trends, and their correlation coeffi- 
cient is ?v(2-14) = 0.55 at the 95% level of confi- 
dence. In contrast, for the recent cycles of n = 15- 
23, there are 3/6 pairs with the same/opposite 
trends, and their correlation coefficient becomes 
negative, at 7y(2-14) = —0.35, at the 63% level 
of confidence. If we consider only the case of 
V min {n) = -1, as for n = 3, 5, 6, 10, 12, 14, 
15, 17, 19, and 23, there are 6/4 pairs with the 
same/opposite trends. For the earlier cycles of 
n < 14, there are 5/1 pairs for V m i n (n) = —1 
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with the same/opposite trends, while for the re- 
cent cycles of n > 15, there are 1/3 pairs for 
V m in(n) — —1 with the same/opposite trends. 
This implies that a decrease in i? m i n is indeed fol- 
lowed by a decrease in R m in most of the earlier 
cycles, while in the recent cycles, a decrease of 
Rmin tends to be followed by an increase of R m . 

In summary, in terms of the varying trend, 
there is no statistically significant positive correla- 
tion between V m and Knin for all cycles (r = 0.18). 
The positive correlation between Vm and V mm ex- 
ists only in the earlier cycles (rv(2-14) = 0.55) at 
the 95% level of confidence. Concerning in the 
recent cycles, this correlation becomes negative 
(rv(2-14) = —0.35). The behavior of the vary- 
ing trend changed in recent cycles. Therefore, we 
cannot infer a decrease of R m from a decrease of 
i? min for Cycle 24. 

4. Temporal Variation in the Running 
Correlation 

In the previous two sections, one has noted that 
the correlation between i? m and R m in behaves dif- 
ferently for different periods of time. Now, we an- 
alyze the temporal variation in the running cor- 
relation with a moving time window of w = 5 
cycles. For each cycle n, we calculate the cor- 
relation coefficient betw een f? m (i) and R m in(i) for 
i = n-2, n-1, . . . , n+2 (|Du et al.ll2009al ). denoted 
by r(5, n). The results are shown in Fig. 2a. 

It can be seen that the correlation is positive 
before n = 13 and significant at the 90% level of 
confidence for cycles n = 5-9 (stars). This implies 
that a lower (higher) level of i? m in tends to be fol- 
lowed by a weaker (stronger) R m for these earlier 
cycles. However, the correlation decreases since 
n = 14, and becomes negative since n = 18, im- 
plying that a lower i? m i n corresponds to a stronger 
i? m (see also n — 15, 17 and 19 in Fig. la). There- 
fore, a lower i? m i„ has not always been followed by 
a weaker R m . In other words, we cannot infer a 
very weak R m of Cycle 24 from the preceding very 
low level of i? m in- 

5. Discussions and Conclusions 

It is well known that the maximum amplitude 
(Rm) of a solar cycle is positively correlated with 
the preceding minimum (R m - m ), so that a low i? m i n 
tends to be followed by a weak R m . However, this 
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Fig. 2. — (a) Running correlation coefficient between 
i? m and i?min for a 5-cycle moving window. Stars indi- 
cate that the relevant values are significant at the 90% 
level of confidence, (b) R m (solid) versa the fitted value 
(dotted) from Equation ||3J). (c) Running correlation co- 
efficient of i? m with both aa m i n and i? m i n for a 5-cycle 
moving window. 



relatio nship is not always effe ctive for individual 
cycles ( Wang fc Sheelev 20091 ). especially for the 
recent cycles, as shown in Fig. 2a. The correlation 
between R m and R mm varies with time (n). 

We analyzed the temporal behavior of this cor- 
relation and the varying trends (V) of i? m and 
i? m ; n . In the recent cycles, they all show a neg- 
ative correlation. Since the prediction of R m re- 
lies more o n the recent cycle rather than on the 
past cycles dSchatten 2005 ; Svalgaard et al. 2005; 
Du et al.ll2008l 2009b ). the negative correlation in 
the recent cycles cannot infer a very weak i? m from 
a very low i? m i n . 

One may argue that i? m and i? m i n have a similar 
shape in the most recent four cycles of n = 20- 
23. Along with the developing trend of these 
cycles, i? m (24) should be very small. However, 
whether this behavior holds true is questionable 
before and after these cycles. It should be noted 
in Fig. la that R m has never decreased in three 
successive cycles. The i? m value decreased two 
cycles from n = 3 to 5, and then leveled off to 
n = 6, and decreased two cycles from n = 8 to 
10, and then increased to n = 11. Now that the 
R m value decreased two cycles from n = 21 to 
23, it seems to increase or level off according to 
its past behavior. On the other hand, i? m i n (24) 
is not the lowest one ever see n. It i s higher than 
cycles 6, 7 (Fig. la), and 15 (Li 2009): i? mi „(24) > 
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(15) > i? mi „(7) > i? min (6). However, cor- 
responding to these local minima, the following 
R m values are not local minima: i? m (6) ~ i? m (5), 
i? m (7) > i? m (6), and R m (15) > R m (U). From 
this information, we cannot yet infer that Cycle 
24 is a local minimum. To say the least, it is un- 
likely that Cycle 24 will be the weakest cycle. 

In conclusion, we have not found sufficient ev- 
idence for the low(est) level of Solar Cycle 24 in- 
ferred from the low level of the present state. The 
sunspot number is highly correlated with other so- 
lar activity indices, such as sunspot group number, 
sunspot area, solar radio flux, and so on. There- 
fore, the above conclusions can also be reached 
when using these indices. 

Near the time of the solar cycle minimum, geo- 
magnetic activity is a much better indicator of the 
ensuin g maximum amplitude (R m ) for the sunspot 
cycle ( Ohll 19661) than the minimum am plitude 
is (Rmxa). lHathawav et al.l (|l999t l2009h tested 
the predictive powers of several methods for cy- 
cles 19-23, and concluded that the geomagnetic- 
related precursor methods outperform the others. 
The minimum smoothed monthly mean aa index 
(aa m i n ) near the time of the solar cycle minimum 
is shown in Table 1, in which the values of cy- 
cles 9-11 are taken fr om the equivalent annual val- 
ues ( Du et al . 2009b). One can note that the vary- 
ing trend (V) of R m follows well with that of aa m ; n 
— with only the two exceptions of cycles 16 and 
22. The correlation coefficient between R m an d 
aa m i„ is usually as high as 0.9 (|Du et al . 2009b). 
The application of aa mm in the prediction of 
can be found , for ex a mple, in Hathaway _ j2009l ) 
and iDu et all (|2009bh . IWilson et al.l (|l998l ) sug- 
gested the bivariate of both aa m - la and i? m i n to 
predict R m . Using the data for cycles 9-23 in Ta- 
ble 1, the bivariate-fit regression equation of R m 
versus both aa min and R m is: 



Rn 



5.0 + 10.56aa n 



3.18i? n 



cr = 15.3, 

(3) 

where a is the standard deviation of the equation. 
Figure 2b shows the observed R m (solid) and the 
fitted R m (dotted) from the above equation. Sub- 
stituting the values of aa m i n (8.4) and i? m i n (1.7) 
into this equation, the peak of the next cycle is 
predicted as R m (24) = 88.3 ± 15.3 (labeled by a 
star). This prediction is close to that predicted by 
the single variate of i? m in in Equation (UJ). But, 



the correlation for the bivariate of both aa m ; n and 
Rm{r = 0.92) is much higher than that of the sin- 
gle variate of i? m i n (r = 0.56). If this prediction 
comes true, Cycle 24 will be modest rather than 
the lowest one. 

The prediction of R m is related to the behav- 
i or of s olar activity in the past cycles 
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(|2009bl ) pointed out that Ohl's precursor method 
performed well only if the related correlation co- 
efficient becomes stronger. If the correlation co- 
efficient becomes weaker, its prediction would be 
questionable. Figure 2c shows the running correla- 
tion coefficient r(5, n) of R m with both aa m ; n and 
Rmin for a five-cycle moving window. It is seen 
that the last value (n = 21, corresponding to the 
data for cycles 19-23) drops drastically. There- 
fore, other methods are needed to check the above 
prediction. 

Predicting the future level of a solar cycle 
is a com plex project in so lar physics and space 



weather (Wa ng et al 



2009). This paper stresses 
that the low level of i? m i n in the present state is 
insufficient to infer a lo w(est) level for Solar Cycle 
24, as suggested byQ ( 2009 ). Whether a predic- 
tion from a simple parameter succeeds is related 
to the behavior of solar activity in the past few cy- 
cles. When a solar cycle is well underway (two to 
three years after the minimum), its behavior can 
be predicte d to a good ext ent with curve fitting 
techniques ( Hathaway 20091 ). 
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Table 1: Cycle maximum (i? m ), minimum (i? m in) and minimum aa index 
(aa m i n ) and their trends (V). 



n 


Rm{V) 


Rmin(V) 




n 


Rm(V) 






1 


86.5 


8.4 




13 


87.9(+) 


5.0(+) 


10.6(+) 


2 


115.8(+) 


11.2(+) 




14 


64.2(-) 


2.7(-) 


5.9(-) 


3 


158.5(+) 


7.2(-) 




15 


105.4(+) 


1.5(-) 


8.2(+) 


4 


141. 2(-) 


9.5(+) 




16 


78.1(-) 


5.6(+) 


9.4(+) 


5 


49.2(-) 


3.2(-) 




17 


119.2(+) 


3.5(-) 


13.2(+) 


6 


48.7(-) 


O.O(-) 




18 


151.8(+) 


7.7(+) 


16.3(+) 


7 


71.5(+) 


0.1(+) 




19 


201.3(+) 


3.4(-) 


16.9(+) 


8 


146.9(+) 


7.3(+) 




20 


110.6(-) 


9.6(+) 


13.8(-) 


9 


131.9(-) 


10.6(+) 


14.1 


21 


164.5(+) 


12.2(+) 


17.2(+) 


10 


98.0(-) 


3.2(-) 


10.3(-) 


22 


158.5(-) 


12.3(+) 


17.5(+) 


11 


140.3(+) 


5.2(+) 


16.0(+) 


23 


120.8(-) 


8.0(-) 


15.9(-) 


12 


74.6(-) 


2.2(-) 


6.7(-) 


24 


? (?) 


1.7(-) 


8.4(-) 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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